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Abstract 

The Casimir effect in an inhomogeneous dielectric is investigated using Lifshitz's theory of electromagnetic 
vacuum energy. A permittivity function that depends continuously on one Cartesian coordinate is chosen, 
bounded on each side by homogeneous dielectrics. The result for the Casimir stress is infinite everywhere 
inside the inhomogeneous region, a divergence that does not occur for piece-wise homogeneous dielectrics 
with planar boundaries. A Casimir force per unit volume can be extracted from the infinite stress but it 
diverges on the boundaries between the inhomogeneous medium and the homogeneous dielectrics. An alter- 
native regularization of the vacuum stress is considered that removes the contribution of the inhomogeneity 
over small distances, where macroscopic electromagnetism is invalid. The alternative regularization yields 
a finite Casimir stress inside the inhomogeneous region, but the stress and force per unit volume diverge 
on the boundaries with the homogeneous dielectrics. The case of inhomogeneous dielectrics with planar 
boundaries thus falls outside the current understanding of the Casimir effect. 
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1. Introduction 

Quantization of the electromagnetic field leads to a vacuum, zero-point energy-momentum that di- 
verges [I]. As this infinite vacuum energy clearly does not exist it was initially assumed that the infinity 
should be zero, until Casimir showed [2] that a finite vacuum energy can sometimes be extracted, leading 
to the prediction of measurable forces. 

In recent years there has been considerable progress in the experimental demonstration of Casimir forces 
between objects [S], and a growing appreciation of the problems and opportunities these forces entail for 
micro- and nano-engineering [4, 5. . On the theory side, it is understood how to calculate the Casimir forces 
between an arbitrary number of separated bodies, the most versatile approach being that of Lifshitz j6j[7l[l] 
wherein the electromagnetic stress tensor and energy density are obtained from the Green tensor for the 
electric field (or vector potential). Although the calculations for any arrangement of objects beyond the 
classic case of two parallel half-spaces [2, [B] are extraordinarily cumbersome, the problem can be solved 
purely numerically using standard codes for the Green tensor [9j 1101 lllj . 

Yet the theoretical understanding of the Casimir effect remains incomplete in some fundamental respects. 
Most importantly, the Casimir self-energy of an object is a quantity whose exact meaning is still debated 
and for which the standard approaches like Lifshitz theory give a diverging answer in general [5]. The 
subject of Casimir self-energy began in 1968 with Boyer's spectacular conclusion that the self-energy of 
an infinitesimally thin, perfectly conducting spherical shell is positive, giving a repulsive (i.e. outwardly 
directed) Casimir force on the shell [12] . An unsettling feature of the spherical shell is the occurrence of an 
additional divergence which is not present in the case of perfectly conducting parallel plates [5] and which 
must be regularized to obtain a finite answer. It might be hoped that the extra divergence is due to the 
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idealized properties of vanishing thickness and perfect conduction of the shell, but the apparently reasonable 
case of spherically symmetric materials with general electric permittivity and magnetic permeability is even 
more ill-behaved. Even with allowance for (temporal) dispersion, the additional divergences remain and 
there is no agreement on their precise significance, or how to remove them p!51 [151 [T51 [T71 [151 [T51 15], 
Hence the striking fact that the seemingly innocuous issue of the Casimir self-energy of a dielectric ball is 
still not understood |15l [3] . Similar problems occur for cylindrically symmetric media, the other case where 
Casimir self-forces have been investigated in some detail. (See [S] for a guide to the Casimir literature on 
spherical and cylindrical geometries.) This lack of understanding of the self- force for these geometries makes 
any attempt to compute the Casimir forces on concentric spherical and cylindrical bodies problematic. 

The work described above has all been for piece-wise homogeneous media, with the non-zero Casimir 
forces arising from the boundaries between the homogeneous components. Although the general theory of 
Lifshitz and his co-workers ^|7J[S] is formulated for inhomogeneous, dispersive dielectrics, to our knowledge 
there is no example in the literature where the Casimir effect for an inhomogeneous dielectric is actually 
calculated|^ In this paper we provide such an example by using Lifshitz theory to calculate the vacuum elec- 
tromagnetic stress tensor for a simple model of an inhomogeneous dielectric. The stress tensor specifies the 
local Casimir self-force in the medium arising from its inhomogeneity. We choose the simplest example where 
the inhomogeneity is in one spatial dimension only, the medium being homogeneous in planes orthogonal 
to this direction. Boundaries are imposed on both sides of the inhomogeneous medium, beyond which the 
material is homogeneous. As the boundaries are planar, we expected to avoid the problems with divergences 
described above, which appear to be associated with curved boundaries and a high degree of symmetry. 
Our results did not meet these expectations. Specifically, we find that the regularization prescription for 
inhomogeneous dielectrics advocated in standard Lifshitz theory |71 Uj yields an infinite Casimir stress ev- 
erywhere inside the inhomogeneous region. A Casimir force per unit volume, given by the spatial derivative 
of the stress, can be extracted; this appears to be finite inside the inhomogeneous region but it diverges on 
the boundaries with the homogeneous dielectrics. We investigate these divergences by considering a new 
regularization procedure, based on a consideration of the limited validity of macroscopic electromagnetism. 
The new regularization gives a finite Casimir stress inside the inhomogeneous medium, however the stress 
and the force per unit volume diverge on the boundaries with the homogeneous regions. We thus find that, 
in addition to spherically and cylindrically symmetric dielectrics with boundaries, inhomogeneous dielectrics 
with planar boundaries contain divergences that fall outside current understanding of the Casimir effect. 

In Section |2] we review Lifshitz theory for a medium inhomogeneous in one spatial dimension. The 
problem of calculating the Casimir stress reduces to solving the equations for two scalar Green functions, 
associated with the TE (transverse electric) and TM (transverse magnetic) modes. We discuss the problem 
of regularizing the quantum- vacuum energy-momentum to obtain a finite Casimir stress. In addition to the 
standard regularization we consider an alternative that removes the contribution to the Casimir stress of the 
inhomogeneity over small distances where macroscopic electromagnetism is invalid. In Section[3]we consider 
a particular functional form for the permittivity of an inhomogeneous dielectric bounded by homogeneous 
regions and obtain analytic solutions for the Green functions. The resulting Casimir stress is considered 
using both the standard and the alternative regularizations; we find that neither gives a Casimir stress and 
force per unit volume that are everywhere finite. 

2. Lifshitz theory for media inhomogeneous in one Cartesian coordinate 

In the formalism of Lifshitz [BJ [7J the forces on macroscopic bodies due to the electromagnetic quantum 
vacuum are obtained by calculating the vacuum expectation value of the Maxwell stress tensor: 



previous study |21) considered layers of homogeneous slabs as an approximation to an inhomogeneous dielectric and 
computed the Casimir energy arising from the TM evanescent modes of the system. 
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The materials are in general inhomogeneous and dispersive with dielectric permittivity e(r,uj) and magnetic 
permeabilty fi(r,uj); anisotropy may also be included but here we consider only scalar s and /i The vacuum 
correlation functions of electric and magnetic fields in (jlj are given by 

/? f°° 

(D(r,t)®E(r,t)) = £(r, ^'G(r, r, z^), (2) 

Jo 

(H(r,t)(»B(r,t)) = - lim /" \ V x G(r, r', i^) x V' . (3) 

In (|2|-(|3|, G(r,r',iO is the imaginary-frequency retarded Green tensor for the vector potential in a gauge 
in which the electric potential is set to zero; its equation is 

Vx ^i— Vx+^£(r,^0)G(r,r^^0 = l<5(r-r'). (4) 

The notation x V' in (3 1 denotes a curl on the second index of G(r,r',i^), so that for a vector V(r') we 

have Vx V'= V x V. Physically, (4) represents a dipole oscillating with frequency i£, at the point r' and 
G(r, r', i^) is the resulting vector potential at the point r. The second index in dj represents the orientation 
of the dipole at r', while the first index represents the components of the vector potential at r. An important 
property of the Green tensor, which holds for systems invariant under time-reversal ^ , is 

G,,(r,r',zO = G,.(r',r,*0. (5) 

The Casimir force F on a body occupying a volume fl is obtained by integrating the stress tensor ([T]) over 
the boundary surface dH. of the body: 

F = [ (T dS = [ V adV. (6) 

Jdn Jn 

An infinitesimal element dV of the body at position r experiences a local force f dV, where f is the force 
per unit volume at r: 

f = V • cr. (7) 

The stress tensor ([T]) is thus sufficient to determine the Casimir forces on the materials, but one can 
also calculate the quantum- vacuum energy density. Since the materials are dispersive, this energy density 
is given by Brillouin's formula |22j . In terms of the monochromatic (in imaginary frequency) expectation 
values 

(E(r)®E(r)>e = -^e'G(r,r,za, (8) 



TT 

(H(r)®H(r))5 = — lim ... V x G(r,r',z^)x V, (9) 



the Casimir energy density is 



. ^ i f (.„«f^^» (E . E), + (fl . AU) . (10) 

The Casimir stress Q and energy density ( 10 1 are thus obtained by solving Q for the Green tensor, 
subject to the boundary conditions imposed by the specific materials. A drawback to this approach is that 
the number of cases where Q can be solved analytically with non-trivial boundary conditions is severely 
limited. On the other hand, the task of solving Q for a given arrangement of objects is a well-posed 
classical electromagnetism problem. It is the issue of regularizing the Green tensor so that it yields a finite 
Casmir force that sets the current limit to standard Lifshitz theory, since it is not understood how to do this 
in general. As described in the Introduction, the Casimir effect for general dielectrics with spherical and 
cylindrical symmetry is still unsolved, even though the Green tensor in each case is known. Before discussing 
the regularization problem for inhomogeneous media, we show how to solve Q in an inhomogeneous medium 
where e(r,(jj) and /x(r,w) depend only on one Cartesian coordinate |531I23]. 

3 



2.1. The Green tensor 

Let £(r,(jj) = e{x,uj) and fj,{r,u!) = fj,{x,uj) depend only on the ^-coordinate. The resulting homogeneity 
in the j/z-plane allows us to Fourier transform the Green tensor: 



G{x,x' ,u,v,i^) — / dy dzG(r,r',i^) 



'iu(y-y')-iv{z-z') 



(11) 



SO that the radiation emitted by the dipole is decomposed into waves with wave- vector components ky — u, 
kz = V at the imaginary frequency i^. With these symmetry assumptions, we are restricted to considering 
planar material boundaries, given by values of x. This is one of the cases where the vector-potential waves 
(proportional to the electric field) can be decomposed into independent TE (electric field in yz-plane) and 
TM (magnetic field in yz-plane) waves. The general solution for the Fourier-transformed Green tensor (111 
can then be written in terms of scalar Green functions gEix,x' ,u,v,i^) and gpiix, x' ,u,v,iS,) for the TE 
and TM waves pSlEij . 

To write the equation for the Fourier-transformed Green tensor (111 we introduce the vector operators 

(12) 



X> and T>', which function as Fourier transforms of the operators V and V': 

f := {dx, iu, iv), T)' (9^, — m, —iv). 
We then obtain from Q and (111 

(-D X — X> X +^e{x,iS,) \ G{x,x' ,u,v,iS) = t5{x~x'), 

and the solution of this equation has the reciprocity property that follows from ([5]) and (111: 

Gij{x, x' , u, V, i^) = Gji{x' , X, —u, —V, i^). 
The scalar Green functions and gM are defined by 

1 - e 



V 



X>- ^e{x,i^) ] gEix,x',u,v,i^) ^S{x-x'), 



I c^^- ^Kx,-iO ] 9M{x,x',u,v,i^) = 5{x~x'), 



and exhibit the same reciprocity symmetry (14 1: 

gE{x, x', u, V, «0 = 9e{x', X, -u, -v, i^), 
gnix, x', u, V, i^) = gM{x\ x, -u, -v, i£) 

We also define a unit vector in the direction of the TE electric field: 



1 




(13) 

(14) 

(15) 
(16) 



(17) 
(18) 



(19) 



The general solution of (13) for the Fourier-transformed Green tensor G can now be written as 

^2 



G 



-ns gE®nE + c 



, T> X ngffM ® ngx T>' ^ 

e{x,iS)e{x' ^iCjS,"^ e{x',i^)C 



;16{x-x') 



s{x',tO^'- 



;nE(E)nES{x-x'). (20) 



To confirm (20 1 as a solution of (13 1 requires the following identities for an arbitrary function f{x) and 
vector V(a;): 



[D X f{x)V X V(a;)]^ = V, f{x)V,V,{x) - V , f {x)V ,V,{x) . 
V-nEf{x) = 0. 
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(21) 
(22) 



The effect of the differential operator in (f3l on the terms in (20 1 containing the scalar Green functions is 
found from the following, which are obtained by successively applying (21), (22 1 and the relevant Green- 
function equation (15 1 or (16): 

V X 



1 ^ ^2 

— — — 1> X 5e = — -e{x, i^) ue qe ~ 5{x - x'), 



[i{x,ii) 



e{x,i^) 



- Ue qm 



1 



fi{x,i^) 



THE S{x — x) 



(23) 



(24) 



When (20) is inserted into the left-hand side of (13) and (23l-( p4|) are employed, one obtains only terms 
containing 5{x — x'); the latter simplify to the right-hand side of (13) because of the identity 



X> X 



fiix.i^) 



■D X [16{x - x') - (g) riij S{x - x')] + X> 



1 



fiix,i^) 



Ue (X" ueSIx - x')x X>'= 0, 



(25) 



which must be verified using (21 1 and the definitions (12) and (19). 

The problem of finding the Green tensor is thus reduced to that of solving the two ordinary differential 
equations ( 15 1 and ( 16 1 for the scalar Green functions. The conditions that qe and 5m must obey at material 
boundaries, values of x where £, fj, or their derivatives are discontinuous, are obtained from (15 1 and (16) 
in the same manner that the boundary conditions for electromagnetic fields are obtained from Maxwell's 
equations |25j . The boundary conditions are continuity of 



1 



9E, 



9M, 



fi{x,i^) 



9x9 E, 



and 



1 



dxgM- 



(26) 



The continuity of ( 26 1 ensures the correct boundary conditions on the Green tensor ( 20 1 , which are familiar 
from basic electromagnetism [5S] since G(r,r',i^), eG(r,r',i^), V x G(r,r',i^), and /i^^V x G(r,r',i^) are 
proportional to E(r,i^), D(r,i^), B(r,i^) and H(r,i^), respectively. 



Written out explicitly using (12), the equations ([15|)-(16l for the scalar Green functions are 

,,2 



1 



rdx 



fi{x,i^) M(a;,«0 
1 



iS,) gE{x, x', u, V, iS,) = S{x - x'), 



-dx 



^ ^fi{x,i^)] gM{x,x',u,v,i^) = 5{x- x'). 



e{x,i£^) e{x,i^) 



In a homogeneous region, where e and fi in (27)-(28l are independent of x, the general solutions are 



ghE = 90hE + ChEie^^ + ChE2e '^^ , 



— A / 'J/2 



y2 _|_ ^2 £(j^)^(i^) 



(27) 
(28) 

(29) 
(30) 



where the cs are arbitrary constants and .gohs, <7ohA/ are the (retarded) Green functions for a space-filling 
homogeneous medium: 



9QhE — 



^ K^O ^~w\x-x'\ 



2w 



gOhM 



— ^(^0 -w\x-x'\ 



2w 



(31) 



In the next section we will solve (|27|-([28|) analytically for an inhomogeneous medium. First we must 
consider the problem of regularizing the Green functions, since use of the full Green tensor in ^ and (|3| 
always yields an infinite Casimir stress ([T]). 
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2.2. Regularization and Casimir stress 

In the formalism developed by Lifshitz and coworkers [51[71[S] for realistic dielectrics, the only meaningful 
quantum-vacuum energy is viewed as that which arises from material inhomogeneity, through smooth or 
sharp spatial variations of e and ^. Thus, according to this viewpoint, the (infinite) quantum- vacuum energy- 
momentum of a space-filling homogeneous medium that emerges from the theory must be regularized to zero. 
This infinite energy-momentum is obtained by using the homogeneous Green functions (31 ) to calculate the 
Casimir stress and energy density (10); hence, at the level of the Green functions, the regularization 
prescription is to discard the homogeneous Green functions (31). 

The simplest case where a meaningful non-zero Casimir energy-momentum arises is that of piece-wise 
homogenous media with sharp boundaries. Here the regularization prescription is clear from that for a 
space- filling homogeneous medium. This is because, for piece- wise homogeneous media, the solution for the 
Green tensor ^ in each homogeneous sub-medium is a sum of the bare or bulk Green tensor, which is the 
solution in the absence of a boundary to the sub-medium, and the scattered Green tensor, which ensures that 
the electromagnetic boundary conditions are satisfied by the total Green tensor on the boundary of the sub- 
medium. (For example, this decomposition into bulk and scattered parts is clear in the scalar Green functions 
(29 1). From what has been said regarding a space-filling homogeneous medium, it is clear that the bulk 
Green tensor must always be dropped, since it has been ruled out as a source of Casimir energy. It follows 
that the Casimir energy-momentum is given by the scattered Green tensor. Unfortunately, the scattered 
Green tensor does not always yield finite Casimir forces on the materials. The most studied examples where 
the scattered Green tensor gives an infinite Casimir force are those of spherically and cylindrically symmetric 
dielectrics. This divergence problem emerged in Boyer's original treatment of an infinitely thin, perfectly 
conducting shell |12j . where he used a mode summation technique. Within standard Lifshitz theory, any 
divergence is supposed to arise from the bulk Green tensor only, and there is no prescription for dealing with 
additional divergences that stem from the scattered Green tensor. Although in a few cases such as Boyer's 
shell, additional regularizations have been performed to obtain a finite Casimir force, there has been no 
proposal for a regularization procedure that will work for general dielectrics with spherical and cylindrical 
symmetry. Moreover, there is no consensus on the exact meaning of the additional divergences or on the 
vahdity of regularizations Uke Boyer's [nj[Tl[Til[Tni[I71[TSl[Il[Mll3]- Methods other than Lifshitz theory 
have also been unsuccessful in obtaining a (finite) prediction for the Casimir force on general spherically 
and cylindrically symmetric dielectrics [5]. 

The paragon of a successful calculation of the Casimir force for realistic materials is the Lifshitz-thcory 
formula for the force between parallel half-spaces made of homogeneous dielectrics, separated by a third 
homogeneous dielectric [7]. This formula has been tested experimentally, both directly and through its 
prediction of the approximate force when one of the half-spaces is replaced by a sphere Here the 
regularization prescription of dropping the bulk Green tensor yields a finite Casimir force on the half- 
spaces, but even in this example the issue of divergences is not entirely straightforward. If the x axis is 
perpendicular to the planar boundaries of the half-spaces, then ([6| shows that the perpendicular force on 
the half-spaces is determined by axx, which is zero inside the half-spaces but jumps to a constant finite value 
between them. Although a^x is sufficient to find the Casimir force, additional information on the vacuum 
energy-momentum is provided by the other components of the (diagonal) cr and by the energy density ( 10 1. 
It turns out that ayy, a^z and p all diverge on the surfaces of the half-spaces. The fact that p is not finite 
seems to raise the possibility of a contradiction in the theory, since the force can also be computed by 
integrating p to obtain the Casimir energy and differentiating this with respect to the distance between the 
half-spaces. If this path to the Casimir force is followed it turns out that the terms causing the divergence of 
p on the boundaries disappear upon differentiation with respect to the distance between the half-spaces, and 
the resulting (finite) Casimir force is exactly that obtained from Uxx ■ The divergence of the Casimir energy 
density on the boundaries of the half-spaces does not affect the Casimir force because it is not caused by 
the cavity between them. A single planar slab of dielectric suffers from the same divergence of its Casimir 



^The following results are probably well known to the experts, but they are not usually recorded in presentations of the 
Lifshitz theory for parallel half-spaces. 

6 



surface energy, a phenomenon noted by Barton, using a completely different theoretical framework, for the 
case of a plasma sheet [26j [27] . Although the Casimir forces on planar homogeneous dielectrics come out as 
finite, the divergences in the energy-momentum tensor cannot be dismissed as unobservable since the total 
energy-momentum tensor has a local significance as the source of the gravitational field. There is therefore 
still an issue of divergences to be explained here, though luckily it does not interfere with the prediction of 
a finite Casimir force. 

Let us turn to the case of general inhomogeneous media, where e and /i are allowed to vary continuously 
in space. In these circumstances the Green tensor no longer has a clear decomposition into a bulk part and 
a scattered part. The regularization procedure advocated in standard Lifshitz theory |7][H] is the following: 
before using the Green tensor to calculate the stress tensor at r, subtract the bulk Green tensor for a 
homogeneous medium whose constant values of e and are equal to e(r) and /i(r). This recipe is designed 
to remove the contribution to the vacuum stress that does not arise from material inhomogeneity. To 
our knowledge, this regularization has never been implemented in a solved example, and in using it for the 
solution presented in the next section we will find that it gives an infinite Casimir stress in the inhomogeneous 
dielectric. The cause of the difficulty may be that, while this regularization excludes the contribution to the 
stress from each infinitesimal (and therefore homogeneous) volume element of the dielectric, it includes the 
contribution made by the inhomogeneity down to arbitrarily small scales where the concept of a dielectric 
permittivity is invalid. Realistic models of e{v^i^) and /i(r,z^) should not vary significantly over distances 
comparable to the molecular size, but if continuous functions of r are chosen there will be some variation of 
e(r, i^) and /^(r, i^) over arbitrarily small distances, and this small-scale inhomogeneity will contribute to the 
Casimir stress if the standard regularization is used. One might therefore seek a regularization that removes 
the contribution to the stress from the small-scale inhomogeneity. Such a regularization can be performed 
by constructing a tensor Go(r, r', that takes some accurate account of the small-scale inhomogeneity and 
reduces to the bulk Green tensor for a homogeneous medium when r' r; subtracting this from the exact 
Green tensor will exclude the small-scale inhomogeneity, hopefully yielding a finite stress tensor. 

To construct the tensor Go(r,r',z^) described above, we specialize to our case of interest, the one- 
dimensional inhomogeneity of the last section. Here we can work with the scalar electric and magnetic 
Green functions gsix, x') and gMix, x'), and we require functions goE{x, x') and ])oMix, x') that take accurate 
account of the small-scale inhomogeneity and reduce to (31 1 when x' ^ x. The medium at small scales is 



probed by waves with large transverse wave-vector components u and w, and we wish to solve accurately 



( 27 1-( 28 ) when the waves reach the scale where the wavelength is much smaller than the distance over which 



e{x,iC) and ^{x,iS) are varying significantly. Looking at (27)-(28l we see that these equations describe the 
regime in question when the following are satisfied: 



d£{x, iS,) 



^{x,i^) dx 



(32) 



w := \ V? + v'^ + £{x, i^)^{x, i^) 



(33) 



Note e{x,iS,) and ^{x,i^) are real for positive imaginary frequencies and for dissipative media they decrease 
monotonically to 1 as ^ ranges from to oo [22J; the condition (32) is therefore unaffected by the 
dependence. Consider (27 1 without the delta function, the solution of which we denote by g: 



—e{x,i£) U = 0, 



(34) 



An approximate solution of (34 1 for slowly varying e{x,iS) and /i(x,i^) can be obtained by the WKB 
method [l^ ; combined with ( 32 1 this yields the approximation 



cj^^ exp(± 



w 



J dx 



(35) 
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Comparing (35) with the homogeneous electric Green function (31 ) and the reciprocity property (17 1, we im- 



mediately obtain the desired goE^x, x'). The magnetic function goAiix, x') follows similarly and we have 



goEKx.x) = exp 

2^ vow' 



~ f \/£{x-,i£.)e{x',i£,) 
goM{x,x ) = „ , , exp 



2/ 



WW' 



' dx 



w dx 



(36) 
(37) 
(38) 



The functions (36)-(37l are a modification of the homogeneous Green functions iSlh that take account of 



the inhomogeneity of the medium over small distances where e(x,i^) and fj,(x,iS,) do not vary significantly. 
This small-scale inhomogeneity is experienced only by waves with very large u and v. We can remove 
the contribution of this inhomogeneity to the Casimir stress by subtracting ( |36p and (37) from the Green 
functions IjE and gM', this gives the regularized Green functions 



ffrogE '■— 9E — goE, 



9m — .9oM ■ 



(39) 



It is clear from the foregoing that the small-scale inhomogeneity can also be removed from the Casimir 
stress by introducing a cutoff for large wave vectors. One might consider doing this by means of spatial 
dispersion, with e and /x going to 1 for wavelengths on the molecular scale; however, wave-vector components 
can only be introduced along homogeneous directions, as we have been able to do in the yz-plane by means 



of (111. A cutoff in the transverse wave vector components u and v will yield a finite Casimir force, but the 



status of such cutoff-dependent results is disputed (see, for example, the previously cited literature on the 
spherical Casimir effect). 

The Casimir force (|6| depends only on the component a^x of the stress tensor, because of the translation 
symmetry in the j/z-plane. In terms of the scalar Green functions, the stress tensor is computed from ([l])-([3]), 
(111, and ( [20I) , with the Green functions in ( 20 1 regularized either by ( 39 ) or according to the standard recipe 
based on ( |31| . Inversion of the Fourier transform ( |TT| ) leads to a double integral du dv, but the 

, so we can set f°° du f°° dv 2n f„°° du 



integrand for axx only contains u and v in the combination u^ 
u^ . The result for ctt 



and 



XX 



IS 



^xx {,x) 



47r2 



00 /'OO 







du / d£,(Txx{x,u,i£^), 







axx{x,u,i^) := u 



2 / 9rcgE , 5rogM \ Iq^- 
W I \ OxOx'grcgE OxOx'gvcgM 

fj, e J ^1 e 



v=0 ' 
x' —X 



(40) 
(41) 



We obtain from ([6]) the (x-directed) Casimir force per unit area Fj^ on a slab of material bounded by x ~ xi 
and X = X2 (xi > xi): 

Fa ^ o-xx{x2) - cTxxixi), (42) 
and the (a;-directed) force per unit volume ([t]) simplifies to 

df^xx 



dx 



(43) 



In the following section we calculate the Casimir stress (40 1-(41 ) for a simple model. As well as employing 
the the standard regularization [8j based on (31), we also show the result when the regularization (39 1 is 
used. 
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Figure 1: The left plot shows the permittivity for a central homogeneous dielectric (e = ec) filling < x < xji, flanked by 
two dielectric half-spaces with e = (left) and e = en (right). The right plot shows the resulting Casimir stress (Txx, zero in 
the half spaces and negative (constant) in between. The Casimir force repels the half-spaces from each other. 



3. Example 

The formalism of the last section was set up for dispersive media, with e{x,i(,) and pi{x,i^) taken as 
functions of (imaginary) frequency ^. As well as being essential for an accurate calculation of the Casimir 
force for specific materials, the inclusion of the measured dispersion properties provides a physical cutofi^ 



in the imaginary- frequency integral in ( 40 1 . We will be encountering divergence problems with the formula 



(40 1, but none of these divergences can be cured by the introduction of (temporal) dispersion since they 
always involve a failure of the integration over the transverse wave vector u to converge. What we require 
is a simple model to test the theory of the Casimir effect for inhomogeneous media; since dispersion does 
not resolve the theoretical issues we do not include it in our results. In addition, we assume the medium is 
a dielectric with no magnetic response {fJ. — 1), so that the inhomogeneous medium is completely described 
by a permittivity function e{x). 

In choosing a function s{x), it is interesting to consider a generalization of an arrangement that gives a 
repulsive Casimir force. The classic configuration of Casimir theory is that of parallel half-spaces made of 
homogeneous dielectrics, separated by a third homogeneous dielectric fT. If the permittivities el, £r and 
Ec of the left, right and central dielectrics obey 



El> £c > £r, 



(44) 



then the Casimir force on the half-spaces is repulsive [7]. (The left and right half-spaces can of course be 
interchanged to the same effect.) The arrangement (44 1 of permittivities is depicted in Fig.[l] together with 
the resulting Casimir stress axx that determines the force. The stress axx is zero in the half-spaces and has a 
constant negative value between them. An integration of the stress tensor over the boundary surface of each 
half space gives the Casimir force ^ acting on it; the non-zero force component Fx is clearly negative for the 
left half-space and positive for right half-space, so they are repelled from each other. The repulsive force per 
unit area on each half space is \axx\- The result is the same if the outer dielectrics have a finite, macroscopic 
thickness, since in the force per unit area (42 1, axx will effectively be zero on the outer boundaries of the 



plates while having the same value between them as for infinite thickness. For real, dispersive dielectrics the 
permittivities will change with frequency and to obtain a repulsive force for a given plate separation, (44 1 
must hold over an appropriate range of imaginary frequencies [7J . (If ( 44 1 holds only for a finite range of ^ 



then, depending on the separation, the force may be repulsive or attractive force when the integration over ^ 
is performed [7]). This remarkable prediction [7] has now been verified experimentally by the measurement 
of a repulsive Casimir force between a gold sphere and a silica plate immersed in bromobenzene |29l . 

Because of the intrinsic interest of the configuration shown in Fig. [l| we consider a generalization of it 
where the permittivity is a continuous function of x so that the central dielectric is inhomogeneous — see 
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X 



Figure 2: A modification of the permittivity in Fig. [TJtfiat cfianges the material in the central region so that e is continuous 
everywhere. With the choice of the displayed exponentially decreasing permittivity in the inhomogcneous region < x' < xji, 
the integrand (UTjl in the Casimir stress ||40| can be found analytically. 



Fig. |2] Our choice of the monotonically decreasing e{x) in the central region was determined by requiring 
that the equations (27l-(28) for the Green functions be solvable analytically. We obtain analytic solutions 



for the Green functions with the permittivity 

e{x) 



She 



— bx 



eL,b positive constants, 



(45) 



and we join this continuously to homogeneous half-spaces on the left at a; = and on the right at x = xr, 
as shown in Fig. [2] This example provides an interesting test of intuition about the Casimir effect, since 
from the result in Fig. [Tjone might guess that the parts of the inhomogcneous medium close to x = and 
X = xji experience a Casimir force that repels them from each other. This would imply that a^x is u-shaped 



in the region < a; < so that the force per unit volume ( 43 1 changes from negative to positive as one 
moves from x — Q to x = x^. One could imagine a solid dielectric with the permittivity of Fig. [2] being 
produced through a manufacturing process that allows the density of the solid to be varied continuously as 
it is constructed; the variation of the permittivity in Fig. [2] would thus be caused by a similar variation in 
the sohd's density. 



We proceed to find the scalar Green functions with the permittivity (451 and /i = 1. As seen from (41 1 



we can set w = 0, whereupon (27 1 and (28) for the permittivity (451 are 



dx^ 



f2 



gs = S{x - x'), 



d 1 



dx El dx 



gM - + 9m = 5{x - x'). 



(46) 
(47) 



Both (46 1 and (47 1 can be recast as the modified Bessel equation with an added delta-function term. Taking 



first the electric Green function ( 46 1 , we make the variable change 

-bx/2 ,^ 



s{x) 



be 



x{s) 



In 



bcs 



The delta function 6{x — x') in (46 1 is related to 5{s — s') by 

5{x{s)~x'{s')) ^ 



\A.[x{s)-x'{s')]\ 



5{s-s') = -bs'S{s-s'), 



(48) 



(49) 
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and the other terms in (46) are easily expressed in terms of s, yielding 



> d'^9E 



dgt 



(50) 



With the delta-function term removed, (50 1 is the modified Bessel equation, independent solutions for which 
are provided by the modified Bessel functions l2u/b{s) and K2u/b{s) ,30J. The general solution of (50) can 
be written 



gs = CElhu/b{s) + CE2K2u/b{s) + ^ [l2u/b{s)K2u/b{s') ~ ^2«/h(s')^2«/b(*)] ~ • (51) 



The term in (51 1 proportional to the Heavyside theta function di^s — s') generates the delta function on the 



right-hand side of ( 50 ) , as is verified by substitution and use of the identities 
d 



^^-/„(s) = ^ [/„-i(s) -I- /„+i(s)] 



ds 



Knis) 



1 



^s/„+i(s) + n/„(s) - ^s/„_i(s) = 0, 
^sKn+i{s) - nKn{s) - ^si^„_i(s) = 0. 



(52) 
(53) 
(54) 



It is important to note that the case s > s', for which the term in (51 ) proportional to 0{s — s') is non-zero, 
corresponds to a; < x', as is seen from the definition (48); similarly, s < s' implies x > x' . After the same 



variable change (48), the magnetic Green function equation (47) can be written 
o d^ (l^ 



s— I ^-gu 
as \ s 



1 



4u' 



1, 



-.9m 



2e 



-S{s-s'). 



(55) 



Without the delta-function term, (55) is the modified Bessel equation for the function gu/s- The general 
solution of ( 55 1 can immediately be written down by comparing it with ( 50 1 and the latter's general solution 



(51) 



_ IkP'ss' 

gM = CMlsIy{s) + CM2SK„{s) + ' [Ii,{s)K^{s') - I;,{s')K^{s)] d{s - s'), 



IF' 



(56) 



(57) 



Equations (51 ) and (56) are the solutions for the Green functions in the inhomogeneous region in Fig. [2] 
The solutions in the homogeneous regions a; < and x > xn are given by ( 29 ) , once we take account that 



the waves from the dipole at x' are outgoing to a; = ±oo; this means that for x < the coefficients of e ""^ 
in the Green functions (29) vanish, whereas for x > x^ the coefficients of e™^ vanish. There are thus four 



unknown constants in the complete solutions for both the electric and magnetic Green functions, and these 
constants are determined by imposing the continuity of ( 26 1 at the interfaces x — Q and x 



xr. 



This is a 



mechanical exercise that yields very lengthy and unenlightening expressions for the constants. The resulting 



exact Green functions are then regularized, either by (39) or by the standard recipe based on (31 ), and the 



integrand (41 1 in the Casimir stress (40) is found. The double integral in (40 1 must be evaluated numerically 
and we proceed to exhibit the results for specific values of the parameters in Fig. [2j 

The first set of parameteter values is depicted in Fig. [3] which shows the permittivity distribution of 
Fig. |2] with El — h — \ and xr = 0.5, so that £r — 1.82. We show first the result of using the 
regularization prescription for inhomogeneous media recommended in standard Lifshitz theory [7,8. Recall 
that the standard regularization is to subtract from the total Green functions, at each x, the Green functions 
for a homogeneous medium whose constant value of e is equal to £{x). Figure |4] shows the behaviour of the 
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1.5 



0.5 



Figure 3: The permittivity of Fig.|2]for the parameter values = 3, b = 1 and xji = 0.5. This gives en = 1.82. 




u=0.5 



20 40 50 100 

u ^ 

Figure 4: The integrand | |41| l in the Casimir stress ( |40[ l with the parameters of Fig.jsjand H = c = 1, where the Green functions 
are regularized according to the standard prescription of Lifshitz theory |3[8|. The left figure plots a^x as a function of u 
with X = O.bxfi and ^ = 0.5, while the right figure plots a^x as a function of ^ for the same x, with u = 0.5. The ^-plot goes 
to zero as ^ — > oo but the n-plot shows a decrease to a constant negative value as -u ^ oo, giving an infinite Casimir stress 
cTxx- Similar plots for other x values inside the inhomogeneous region < a; < xji show the same qualitative behaviour; the 
standard regularization therefore gives an infinite Casimir stress ctxx throughout the inhomogeneous medium. 



integrand axx if this regularization is used. The figure plots the resulting axx as a function of u with ^ fixed, 
and as a function of ^ with u fixed, for x = O.bxR. The u-plot for fixed ^ approaches a constant negative 
value as u — > oo, giving an infinite value for the integral (40 1 at this value of x, and this is the behaviour 
for any x inside the inhomogeneous medium < a: < xr. The standard regularization thus gives an infinite 
Casimir stress axx throughout the inhomogeneous region. 

One might consider falling back on the fact that the Casimir force ( 42 1 on a slab of the dielectric is 
given by subtracting the stress at the two boundaries of the slab, and that a subtraction of two infinite 
values of the stress may give a finite answer. But one is free to consider a slab with one boundary in 
the inhomogeneous region of Fig. [3] a nd the other boundary in one of the homogeneous regions. With the 
standard regularization the force (42 1 on any such slab is infinite, since axx in the homogeneous regions is 
zero. Moreover, the legitimacy of an additional regularization to remove the infinite stress visible in Fig. [4]is 
questionable, since according to the standard theory [71 18] there should be no such infinity. As discussed in 
some detail in the last section, the expedient of additional regularizations on a case by case basis is a matter 
of continuing debate and disagreement. Nevertheless, for information purposes we consider the maneuver of 
taking the derivative in (43 1 inside the integrations in the infinite axx] we can then differentiate the finite 
integrand axx and afterwards perform the integrations to obtain a value for the force per unit area /. This 
could remove an x-independent infinity in axx and give a finite expression for /. Figures [S] and |6] show the 
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Figure 5: The derivative of the integrand with the parameters of Fig. [s] and h = c = 1, where the Green functions are 
regularized according to the standard prescription of Lifshitz theory [71[8]. The left figure plots da^x/dx as a function of u for 
four values of x, with ^ = 0.5, while the right figure plots da^x/dx as a function of g for the same x values, with u = 0.5. It is 
clear that the plots for a: = and x = xn fail to approach zero for large u or for large ^, so integration of daxx/dx with respect 
to u and ^ yields diverging values for the force per unit volume / at these two positions. The values of x chosen in the plots 
are close to or on the interfaces a; = and x = xr; similar plots for x near the center of the inhomogeneous region give curves 
with pronounced oscillations, as shown in Fig. [6] 
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Figure 6: The same as Fig. [s] but for the single position x = 0.5x/j at the centre of the inhomogeneous region. The curves are 
more oscillatory that those in Fig. [5] for positions close to or on the interfaces x = and x = xr. 



behaviour of da^x/dx at different values of x. It is clear from Fig. [5] that da^x/dx at the interfaces x = 
and x = xr = 0.5 is such that its integration would give a diverging value for the force per unit area / at 
these positions. The plots of dcixx/dx inside the inhomogeneous region < a; < xn look more promising in 
that they approach zero for u ^ oo and ^ — > oo. However, we find that the curves for duxx/dx as a function 
of ^ approach zero rather slowly as ^ — > oo. Our attempts to numerically integrate dcfxx/dx for a range of 
a; in < a; < xn produced values that varied significantly with small changes of x and it is not clear that 
the integrals converge. In any case the obvious fact that the integration of duxx/dx gives a diverging force 
per unit area on the interfaces x = and x = Xfj shows that the standard regularization does not give an 
acceptable solution to the problem. 

We now show the results for the alternative regularization ( 39 1 based on ( 36 1 and ( 37 ) . Figures [7]and|8] 
show plots of the integrand axx as a function of u with ^ fixed, and as a function of ^ with u fixed, for 
several values of x. It can be seen how the area under the curves in the u-plots increases significantly as 
X moves towards or xr = 0.5, but the curves approach zero as u ^ cx) and ^ ^ oo for < x < x^, 
giving a finite answer for the stress axx at these values of x. The behaviour of function of u when 

a; = and x — xn, however, shows a divergence of the stress axx at these positions, since axx increases to a 
constant positive value as m — > oo, giving an infinite value for the integral (40 1. Figure [9] shows the resulting 
Casimir stress axx, obtained by numerically integrating (41 1 in (40 1. The integrations must be performed 
separately for each x and the dots show axx for a discrete set of x values, the results for which are joined 
together by a polynomial fit. In the homogeneous material at a; < and a: > 0.5 the stress axx is identically 
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Figure 7: The integrand | |41| l in the Casimir stress | |40[ l with the parameters of Fig. [s] and h = c = 1, where the Green functions 
are regularized by (jsojl and | |36| l- |37| . The left figure plots axx as a function of u for three values of x, with § = 0.5. For 
X = xji = 0.5 the integrand attains a constant positive value as « ^ oo, giving an infinite Casimir stress a^x at this position. 
The right figure plots axx as a function of § for the same three values of x, with u = 1; the curves approach zero as ^ — + oo. 
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Figure 8: The same as Fig. [7] but with values of x less than or equal to 0.5a;^ = 0.25. For x ■■ 
similar to that for x = xr in Fig. [7] giving an infinite Casimir stress axx at this position. 



: the behaviour of the M-plot is 



zero, as we verified by calculating the integrand (41) in these regions where there is no difficulty with the 
standard regularization. Figure [9] also shows the corresponding Casimir force per unit volume (43); in the 
inhomogeneous region / was found by differentiating the polynomial fit to a^x- The unlimited increase in 
a XX as the interfaces a; = and x = xr are approached can be seen in Figure [9] leading to similar divergences 
in /. 

The expectation of a u-shaped axx, giving a Casimir force / that changes direction as one moves across the 
inhomogeneous medium, is to some extent borne out by Fig[9] It was the failure of the standard regularization 
recipe to give a finite stress and force per unit volume that led us to try the alternative proposal based on 
(36) and (37), an approach that at least has some a priori physical justification. Although the alternative 
regularization gives a finite Casimir stress inside the inhomogeneous medium, the stress and the resulting 
force per unit volume are still infinite on the interfaces with the homogeneous regions. We have therefore 
not obtained a physically acceptable solution to the problem considered, with either regularization. 

An interesting parameter value to investigate is er = 1, so that the inhomogeneous medium in Fig. |2] 
joins continuously to vacuum on the right-hand side. An example of this is plotted in Fig. [TO] where the 
other parameters are chosen as El = 3 and b = 1, giving xr = ln(3). The results for the stress axx and force 
per unit volume / in this case are qualitatively similar to those for the previous parameter set. Figure 11 
shows axx and / computed using the alternative regularization based on ( 36 1 and ( 37 ) . The Casimir force 
again diverges on the boundaries between the inhomogeneous medium and the homogeneous regions, even 
though the density of the inhomogeneous dielectric goes to zero at the right-hand boundary. 
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Figure 9: The Casimir stress l |40| (left) and Casimir force per unit volume ( |43| (right) for the permittivity in Fig. [s] with 
h = c = 1, obtained using the regularization given by by | |39| l and | |36| l-| |37l l. The dots in the plot of a^x show the results of 
numerical evaluation of the double integral in \A0\ at those values of x; the line through the dots is a polynomial fit. In the 
homogeneous regions x < and x > 0.5 the stress cr^x is zero. The force per unit volume / in the inhomogeneous medium was 
obtained by differentiating the continuous fit to axx- 
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Figure 10; The permittivity of Fig. |2]for the parameter values = 3, = 1 and 6 = 1. This gives xji = ln(3). In this case 
the homogensous medium on the right-hand side x > xji is vacuum. 

4. Discussion 

We have presented a calculation of the Casimir self-force in an inhomogeneous dielectric, using a simple 
model for the dielectric permittivity. As far as we have been able to ascertain, this is the first analysis of 
the Casimir effect for inhomogeneous media, as opposed to piece-wise homogeneous media. We used the 
formalism of Lifshitz, which is routinely recognized as the most versatile and general theory of Casimir 
forces. Although the standard Lifshitz regularization prescription was formulated with the general case of 
inhomogeneous dielectrics in mind [7, 8 , it gives a infinite Casimir stress everywhere inside the inhomo- 
geneous medium in our example. An attempt to extract a finite Casimir force per unit volume from the 
diverging stress failed. We also tried a new regularization method, which aims to remove the contribution 
to the Casimir force arising from the inhomogeneity over short length scales where the use of macroscopic 
electromagnetism is unphysical. The new regularization gives a finite Casimir stress inside the inhomoge- 
neous medium in the example considered, but the stress and force per unit volume increase without limit 
at the boundaries joining the inhomogeneous dielectric to homogeneous regions. Thus, even with the new 
regularization, we have been unable to solve the example using the accepted theoretical machinery. 

It is probable that if the sharp corners in e at a; = and x = xr in Fig. [2] are smoothed out sufficiently 
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Figure 11: The Casimir stress l |40| (left) and Casimir force per unit volume (jJsjl (right) for the permittivity in Fig. |10| with 
h = c= 1, obtained using the rcgularization given by by ||39[ l and l |36[ |-| [37l l. As in Fig. [9] the dots in the plot of (t^x show the 
results of numerical evaluation of the double integral in HoF. The force per unit volume / in the inhomogeneous medium was 
obtained by differentiating the continuous fit to a^x- 



then the new regularization procedure will give a finite Casimir stress Gxx everywhere. This certainly seems 
to be the case for a permittivity profile given by (45) for all x, although we do not discuss this example here 
since e — s- cx) for x — > —00 and e — > for a; —> 00 is not a very instructive model. In this regard it is important 
to reflect on the fact that the new regularization method achieved a finite stress inside the inhomogeneous 
region by removing the contribution of small-scale inhomogeneity. The sharp corners in Fig. [2] are clearly at 
variance with an approach that views small-scale inhomogeneity in e as unphysical and the regularization 
procedure does not apply at the corners themselves since they are boundaries between regions where the 
stress is regularized separately. A model where the function e{x) is C" (n > I) everywhere and constant for 
a; < and x > Xfi would probably have to be solved entirely numerically. Another case of interest would be 
a function e{x) with discontinuities where de{x)/dx — > at each side of the jump, as in Fig. [l] it is unclear 
if the new regularization would always give a finite stress in these circumstances. 

On the other hand, it is conceivable that the divergence of the regularized Casimir force is not an artifact 
of the theory. The Casimir force could indeed tend to infinity at discontinuities in the derivatives of the 
dielectric profiles, for example at the sharp corners in the graph of the dielectric function of Fig. 2. In 
response to such a strong force, the dielectric may either consolidate a sharp boundary between dielectric 
layers, where the regularized Casimir force is finite, or smooth out the transition region. But of course, these 
ideas are speculations at present. 

As discussed earlier in this paper, it has been known for some time that Lifshitz theory also founders for 
general spherically or cylindrically symmetric dielectrics, and the Casimir force in these cases is unknown. 
A legitimate line of enquiry is whether a Casimir self-force on a sphere or cylinder, or on an inhomogeneous 
dielectric, is measurable. This consideration may play a part in resolving the difficulties, but such materials 
certainly alter the local quantum vacuum and this affects the total stress-energy-momentum of the system. 
The total energy-momentum tensor of a material has a local significance as the source of the gravitational 
field, so it needs to be explained exactly how this quantity is to be computed, for any material, and this 
requires an understanding of the Casimir effect that is currently lacking. Lifshitz theory takes a semi-classical 
approach to the Casimir effect, wherein the electromagnetic field is quantized but the matter is treated 
classically. It may be that the regularization difficulties that have been encountered here and elsewhere are 
due to an inherent limitation of this semi-classical approach. Alternatively, a better understanding of the 
physical basis of the regularization of quantum- vacuum energy may resolve the issue, allowing a prediction 
of the Casimir force in all situations where the required Green functions are known. A final, encouraging 
thought is that experimental investigation of Casimir forces is a maturing discipline ^ , so there is reason to 
hope that experiments can be designed that will shed some light on these areas where the theory is sorely 
underdeveloped. 
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